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SUMMARY 


Motion system drive philosophy end corresponding real-time softv/arc 
have been developed for the purpose of simulating the characteristics of 
a typical synergistic Six-Post Motion System (SPMS) on the Flight Simulator 
for Advanced Aircraft (FSAA) at flASA-An.cs which is a non-synergi Stic motion 
system. This paper gives a brief description of these two types of motion 
systems and the general method of producing motion cues on the FSAA. An 
actuator extension transformation which allov/s the simulation of a typical 
SPMS by appropriate drive washout and variable position limiting is des- 
cribed. 
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1.0 INTRODUCTION 

The use of a moving-base flight simulator to supply motion cues to 
the pilot requires an organized scheme of coordination, washout, and 
limiting of drive coir.mands to produce realistic cues while remaining 
within the physical constraints of the motion system. The motion drive 
software for both the SPMS at NASA-Langley and the FSM at NASA-An:es is 
based on the general techniques developed by Schmidt and Conrad [1 and 2]. 
Hov/ever, the design of the software for these two motion systems is highly 
dependent on the specific system's characteristics. 

A more involved problem is that of developing a scheme for producing 
motion commands to the FSAA resulting in motion cues similar to those 
produced by a typical synergistic SPMS. Such a scheme v/as required for 
the purpose of evaluating different motion systems during a joint NASA/ 
USAF experiment recently conducted at Ames. 

This paper includes a discussion of the characteristics of a typical 
synergistic SPMS and the non-synergistic FSAA. The motion drive software 
for the FSAA is reviewed and the scheme for simulating a typical SPMS on 
the FSAA is presented in detail. 
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2.0 nFSCRIPTION OF KOTION SYSTEMS 


All moving-base aircraft simulators possess position, velocity, and 
acceleration constraints which prevent them from exactly reproducing the 
full motion histories of an aircraft. These constraints vary with the 
number of degrees of freedom and performance limitations of any specific 
motion system design. 

The FSAA at NASA-Ames consists of a transport- type cab integrated 
into a system of drive mechanisms such as electrical motor/generator 
combinations, drive belts, timing chains, and gears which form a six- 
degree-of-freedom motion system [3]. A photograph of this system is 
shown in Figure 1. Although a detailed exposition of the performance 
capability of the FSAA is not presented in this paper, a summary of this 
capability is given in Table 1. The FSAA is capable of motion in each 
axis independent of the other axes (i.e., a non-synergistic system). The 
current FSAA drive technique incorporates a velocity command scheme with 
position, velocity, and acceleration sensors supplying feedback to the 
software or the motion monitoring system. The most unique feature of 
the FSAA is its + 40. feet of lateral travel capability. 

A typical SPMS consists of a cab mounted on a platform v^hich is 
driven by six servomechanisms utilizing hydraulically pov;ered actuators 
or legs. A photograph of such a system (the one at NASA-Langley) is 
shown in Figure 2. A schematic representation of the actuator geometry 
for this system is given in Figure 3. The points where the actuators 
connect to the moving platform form an equilateral triangle and the 
actuators have a minimum length of 8.60 feet (2.62 m) and a maximum 
length of 13.58 feet (4.14 n). When all six actuators are at their 
neutral position (midway between minimum and maximum extension) the 
platform achieves a neutral height (Z*) of 8.33 feet (2.54 m). This 
actuator geometry and the performance capability presented in Table 2 
(not exactly the capability of the NASA-Langley system) were considered 
to be typical for the SPMS model used later in this development. A 




1 


typical SPMS does not have drive systems for each degreee-of-freedom but 
achieves motion in all six degrees by a combination of actuator extensions 
(i.e., a synergistic system). The transformation which converts simulator 
position information to leg extensions for this type of system is discussed 
in detail in this paper (Section 4.0). 
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3 . 0 MOTION SYSTEM DRIVE PHILOSOPHY FOR THE FLIGHT SK^ULATOR FOR ADVANCED 
AIRCRAFT (FSAA) 


Designing the drive logic for any motion system involves establishment 
of a scheme by which motion cues are transmitted to the pilot while keeping 
the movement of the simulator within its constraints. After a motion cue 
has been transmitted, tiie motion system should return to its neutral posi- 
tion without the pilot being aware of this movement. The resultant tendency 
to keep the simulator near its neutral position m.aximizes the allov.able 
motion for subsequent cues. 

The motion system drive philosophy for the FSAA involves producing 
coordinated motion cues by generating high frequency translational and 
rotational accelerations supplemented by low frequency forvjard and side 
-forces via proper cab rotation (i.e., orientation of the pilot within the 
gravitational field). The latter technique is corrimonly referred to as 
"residual tilt". The major conjecture in this philosophy is that very low 
frequency accelerations are not sensed and/or are not important to the 
pilot with respect to motion cue adequacy for most tasks of interest. 

Other factors involved in the drive philosophy are discussed in detail in 
[4]. 


A simplified schematic description of a typical sin.ulation configura- 
tion on the FSAA is shown in Figure 4 and a conceptual block diagram of 
the motion system drive logic is presented in Figure 5. This diagram 
presents the relationship between longitudinal translation and pitch motion 
or the relationship between lateral translation and roll motion, since 
these motion pairs are coupled in the same manner. The vertical transla- 
tion and yaw motion of the simulator are treated independently. There is 
no "residual tilt" term for the yaw axis and no compensation in the 
vertical axis due to high frequency yaw rotations. 

For the translational axes the pilot station accelerations are scaled 
and passed through second order high pass filters. High frequency transla- 
tional acceleration components arc then integrated to obtain washed out 
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translational velocity components. Similarly, rotational velocities are 
■'obtained in roll, pitch, and yaw. The rotational velocities, however, 
are still in the aircraft body axes. Gimble angle transformation from 
aircraft tody axes to simulator cab axes produces simulator drive rota- 
tional velocity commands. These rotational rate commands are modified 
by "residual tilt" terms ccm.puted from low frequency translational 
accelerations. In order to give the effect of coordinated motion cues 
the forward and side forces resulting from washed out high frequency 
rotational rates are subtracted from the high frequency translational 
cues. The final translational rate commands are transformed from air- 
craft body axes to simulator cab axes. Both the translational and 
rotational velocity cotnnands are software limited. To keep the simu- 
lator from drifting, position feedback is used to compensate the out- 
going velocity coinnands to the simulator. 





4.0 ACUTATOR EXTENSION TRANSFORMATION FOR TH E SYNERGISTIC SIX-lOSjjCTOl 
SYSTEM (SPMS) 


To obtain actuator lengths from a given platform centroid position 
a transformation is required. This transformation [5] transforms a given 
vector of three translational positions and three rotational orientations 
into six leg lengths of the six actuators. 

The base platform and the moving platform are shovm in Figure 6. 
Consider each platform of the six-legged motion system to have an associ- 
ated coordinate system as shown in Figure 6. The fixed reference frame 
(x . y , z ) has its origin at the centroid of the fixed platform and the 
moJing°reference frame(x^, y^, zj has its origin at the centroid of the 
moving platform. The two reference frames are both right-handed coordi- 
,nate systems and their axes are aligned when the moving platform is at 
its neutral position. For any actuator (i) the vector relationships 
between the two reference frames are illustrated in Figure 7 and can be 

expressed as 


r, • 

f, = 


( 1 ) 


i = 1 , ... ,6 


which yield! 


£. = A.-B.+R (2) 


A graphical representation of these vectors for all six actuators 
is presented in Figure 8. 


Equation (2) gives the orientation and length of each actuator as 
a function of the dimensions of the fixed and moving (or payload) platform 
and the orientation of the m,oving platform with respect to the fixed plat 
form. Equation (2) is expressed strictly in the fixed platform reference 
frame. Hence, the vector B. is known from the dimensions and leg 
attachm.ent geometry of the fixed platform (Table 3). However, the vectors 
A.j dnd R dre yet to be deterrr.ined. 


If we define (T^, k^) as the triad of unit vectors in the fixed 

platform reference frame, we can represent A^, and R as 


t. = Jc^ i j +£^ k 
1 10 10 10 






B. = 


3 ^+E^ k^ 

10 10 10 


R = R" j\+R^ 


It will be shown later that the components (r’^, 
in terms of simulator coordinates (x^, y^, z^). 


> (3) 

J 

R^, R^) can all be expressed 


Defining (i^, k^^^) as the triad of unit vectors in the moving plat- 

form reference frame and (<», e, as the ordered rotation in yaw, pitch, 
and roll of the moving frame with respect to the neutral position in which 
the moving and fixed frames are aligned (i.e., the moving axes are parallel 
to the fixed axes), the relation betv/een the two reference frames is given 
by 




— 

^0 



o 

1 

= [T] 


L'oJ 


3 

1 


(4) 


where [T] is the Euler angle transformation matrix and is expressed as 


[T] 



^12 ^13 
^22 ^23 
^32 ^33 


(5) 
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and 


= cosv'cose T^2 " cosvsinosinj ' cos^sinecos^; 


-sin;cosi> 


+sin;sin* 


= sin;cose T22 = sini^^sinosin? T22 = sin;sinecoss*> 

+cosi;'Cos<j -cos’i^sinj 


( 6 ) 


T3, . -sino 


1^2 = cosesin<{) 1^2 = cosecos<|) 


The vector from the centroid of the moving platform to the attach 
point of any actuator (i) can be expressed in the moving frame as 


1 ,m 1 ,m m ^ ,m m i ,tn m 


(7) 


where the components (A^ ,aY ,A^ ) are knov/n from the geometry 
of the moving platform (Table 3). The A- vectors in the fixed frame are 
then given by 




( 8 ) 


and equation (2) can now be written as 


'i = [T] 


( 9 ) 


To obtain the scalar form of the actuator extension transformation 
given by equation (9) we multiply 2,^ by its transpose iT giving 

S = N,„-VR)T([T] 

= 'fT] Ai.„,-B/R) 

= A{_JTf[T] A,_„-Si_^[lfD,4Aj_„LT]'R-B{[T] A,_„,4s{3, 

-bTr r''^[Tj A. -R^B.+R^R 
1 1 ,m 1 
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Equation (13) can be simplified to yield 

*(R^)'*2RnA^,„T2,*A^_j22^A^,,„T23-B^) 

+(r" (A^ _„Tj, +A^ _„T 32+A? T33-B^ ) 

-2B1(A’;,Jii^A>,,T,2^A^,J,3) 

-2B1(A^,J2,-A^,J22^A^.J23> 

(A? , J31*a 1 ,j32*Ai , J 33 ' 


( 14 ) 


Equation (14) is still in tcras of i” , r'' and R^. These quantities 
are still unknown and need to be related to the translational position of 
the simulator. In order to do this a transformation between the centroid 
of the moving platform of the motion base at its neutral height and the 
origin of the linear position follow-up measuring potentiometers is neces- 
sary. The coordinate system (i^, k^) to measure the position of the 

simulator for the FSAA is a left handed coordinate system whereas the 
coordinate system (i„, 3„. k„) and (i„, 3„. k„) are right handed coordinate 

systems. 

The different coordinate systems are shown in Figure 9. For any 
arbitrary position of the simulator given by s defined as 


the relationship between the components of the vectors R and S are given 
as 



> (16) 

and = -(z^+z *) J 

z* = neutral height =8.33 feet 


where 


5.0 DRIVE PHILOSOPHY FOR SIMULATING AN $PM$ ON THE FSAA 

The drive philosophy adopted for simulating a typical SPMS on the FSAA 
v/as to preserve the general structure of the FSAA philosophy (Figure 5) as 
much as possible but to incorporate the synergistic effects of an SPMS by 
variable position limiting. This philosophy is somev;hat different than one 
that might be used for designing drive logic for a real SPMS. An example 
of a real SPMS philosophy can be found in Reference 6. 

A conceptual block diagram of the scheme developed for simulating an 
SPMS is presented in Figure 10. 


The current position of the FSAA is given by the vector (x , y , z , 
i>» 0» <{') at any instant of time as measured by the linear and rotational 
position follow-up potentiometers. Conceptually, had there actually been 
a six-actuator system the position vector (x^, y^, z^, ({., e, 4,) would have 
resulted from a certain specific combination of six-actuator lengths. 

Hence, the measured position (x^, y^, z^, a, 0, .s) is used in the actuator 
transformation equation given by equation (14) to determine the six-actuator 
lengths. The actuator transformation equation requires the six-actuator 
motion base dimensions components of the Euler angle transfor- 

mation [T^, and the coordinate transformation between the vectors R and S 
(equation (16)). These computations are done prior to using the actuator 
transformation equation as shown in Figure 10. The s ' actuator lengths 
are then compared with each other and using the longest and the shortest 
actuators a prediction is made of the translational position limits (x^^, 
z^) based on current position of the simulator. For the rotational 
axes, fixed angular excursion limits are used. These six 
limits restrict the subsequent motion cue that can be transmitted to the 
pilot. The details of the method of predicting the position limits are 
given in Appendix A. The desired position commands from the FSAA washout 
scheme are then compared to the six position limits. Based on this com- 
parison the objective is to compute velocity limits. As shown in Figure 10 
if a position limit is exceeded then the velocity limit is set equal to zero. 
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otherwise, a position-velocity parabolic limiting is performed. This 
parabolic limiting scheme is discussed in detail in Appendix B. The 
desired velocity cominand as obtained from the FSAA motion washout is then 
limited to the computed velocity limit. As a final limiting, this limited 
velocity is compared with the past velocity command to predict the result- 
ing acceleration. The final velocity command is then limited to produce 
accelerations within the lower acceleration capability of the two motion 
bases. The outgoing velocity ’s then compensated to prevent the simulator 

from drifting. 


6.0 CONCLUDING REMARKS 


This report has presented a technique of modeling synergistic 
motion simulation on a non-synergistic motion simulator. Such motion 
simulation software can also be used for simulation of different types 
of actuator driven m.otion simulators for evaluation, testing and re- 
search purposes. 

The report has discussed and compared the characteristics of a 
synergistic actuator driven motion base and the FSAA which has a non- 
synergistic independent axis driven system. The general problem of 
the motion drive logic design has been discussed and the washout scheme 
for the FSAA has been outlined. The limitations and the restrictions 
of an actuator driven motion base have been incorporated into the motion 
drive scheme of the FSAA. 

This did not require major modifications to the FSAA washout scheme. 
In particular, an actuator extension transformation is required but no 
iterative inverse transformation is required. 

Integrating this drive scheme within a recent KC-135 aircraft 
simulation required some modification to the washout characteristics in 
order to match the SPMS limiting envelope. However, the general opinion 
was that the simulated SPMS gave the effect of an actuator driven motion 
system. 

It should be noted that only synergistic systems with performance 
capability within that of the FSAA can be simulated with this scheme. 


APPENDIX A 


METHOD FOP PREDICTING SP MS POSITIGM LIMITS 


It has t'-en shov,n that the lengths of the actuators It, I for a given 
position (x,. V . z,. f. 0, 0) of the o«aving platfom portion of a typical 
SPMS with actuator attachn'ent geometry A, and B, can be expressed as 

*4*^^5<''l,m^21*''Un/22*''i,ni^23‘®P 

t(Zj+z*)^-2(Zjn*)(A* ,„T3,iA?_j32tA!_j33-B?) 


-2B^(Alm^21*''i.mT22*''lm^23) 

-2Bl(A^.j3l"''i,nV''i.mT33> 


(Al) 


i ~ l5-*#»6 

These lengths can be used to predict the remaining travel which is crucial 
in the determination of maximum motion cue generation capability of the 
SPMS under consideration. 

The displacement of a typical SPMS in any one degree-of-freedom 
alters the maximum displacement that may then be achieved individually 
in each of the remaining five degrees-of-frecdom. Since there can be 
almost an infinite number of possible combinations of the six displace- 
ments for this type of system, there exist an egually large number of 
position limits. Hence, the set of position limits of an SPMS is not 
fixed and must be dotcriiinad for each set of displacements. 
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The following technique for predicting the position limits of an SPMS 
is similar to the one used at flASA-Langley on their synergistic motion base 
[6]. This technique assumes tfiat the most strict constraint on travel, at 
any instant of time, is governed by the position of the shortest and longest 
actuator at that instant. 

Having obtained the current six actuator lengths using equation (Al), 
the maximum and minimum length actuators can be identified. The current 
position limit for each axis is then determined by fixing the other five 
displacements at their current values and choosing the value of the dis- 
placement under consideration that is the minimum of (1) tiie value speci- 
fied as an absolute limit, (2) the value which would produce maximum exten- 
sion of the longest actuator, or (3) the value which would produce maximum 
contraction of the shortest actuator. 

As an example, consider the prediction of the current position limit 
for x^. From equation (Al) we determine j such that 

- l^il , i = 1,...,6 (A2) 


substituting jljl for \l^\ in equation (Al) we obtain 

|l.|2 = x2.2x^G,.Ej 


(A3) 


where 6j and Ej are dependent on the other five degrees-of-freedom (y^, z^, 
(}>, 0, i|») at that instant and can be expressed as 


G. = A’^ T,,+A^ T,,+A^ 
j j,m 11 j,m 12 j,m 13 j 


(A4) 
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I 


f 


and 


E,- = y^+2y^^^. T.,+A^ T-,,+A^ To^-B^) 
j -^s •'s' j,m 21 j,n 22 j,m 23 j' 






2 ,/aZ \2 z~X\2 , / r,y ./r>Z\2 


-2B^(A^ Ttt+A^ T,o+A^ T,,) 
J J.m n j,ni 12 j,m 13 


-2B^(A^ Joi+A^ T,,+A? T,,) 
J J>m 21 j,m 22 j,m 23^ 


2B^(a’^ T,,+A^ T-.+A^ T,_) 
J J.m 31 j.tn 32 j,ni 33' 


(A5) 


Equation (A3) is a quadratic equation in x . Setting |£.l equal to its 

^ J 

maximum possible value, equation (A3) can be solved to yield the value 
for v/hen the longest actuator has reached its full extension. Because 
of the notion limitations of the base no y^, z^, o, v combination can 
be obtained that gives two roots of equation (A3) of the same sign. The 
positive root is chosen if the x^ velocity is positive. This causes the 
longest actuator to be fully extended with x^ increasing in value toward 
the positive x^ position. Similarly the negative root is chosen if the 
x^ velocity is negative. 

Next, from equation (A1 ) we determine k such that 

1\1 < l^il i = 1,...,6. (A6) 

By setting IL [ to its shortest possible length and solving an equation 
like 

Xg+2V^Ek = (A7) 

where and are similar to and Ej except for k replacing j. The 
quadratic equation (A7) is solved and the solution for the shortest 
actuator is chosen in the same manner as for the longest actuator. 

- 18 - 




✓ 


Hence tv/o predicted values of the maximum available excursion are 
obtained. One is based on the longest actuator reaching maximum extension 
and the other is based on the shortest actuator reaching minimum extension. 
These two values are compared with a specified absolute limit for the 
system and the smallest of the three is chosen as the predicted x^ limit. 
In case of imaginary solutions for the quadratic equations the limit is 
chosen as the specified absolute limit; i.e., imaginary solutions are 
ignored. 

The same method is used for y^ and z^. However, no prediction is 
attempted for angular position limits. These limits are fixeo for the 
system under consideration. 
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APPENDIX B 
POSITION-VELOCITY PAR/\Si 


Consider that any given axis § of i 
given by 

Position limit 
Velocity limit 
Acceleration limit 

Ignoring acceleration constraint, 1 
within 1§1 < and 1§[ < which is il 



Figure P.l: Position-Vel 
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By definition, 


! _ d§ 
dt 


(Velocity) 


^ - dtMt^ 


d d Q5 

dt^§) - 


l-j|- (Acceleration) 


If we wish to establish the conditions for a fixed (or limiting) 


acceleration. 


§ = +a, 


Then equation (B2) yields 


^ d§ 


+a^d5 = §d§ 


integrating. 


+a, d§ = /• § § 


where 


5 = X, 


+a^(x-x^) 




Since we wish to have zero velocity at our position limits, the 
boundary conditions are 


X, = iPt 


X, =0 
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Then equation (B7) yields 


1 


+a^(x+Pj^) ^ j k 


(B9) 


or 


l_ ‘2 

‘i 


X = t^;X ±p, 


(BIO) 


This expression contains the follov.'ing four distinct possibilities 


* - P. - 257 

-x^ 

* - P» ' 2a, 

a 

* P. - 257 
. -;2 
« * P. = 257 


(BID 

(B12) 

(B13) 

(BI4) 


These four equations are represented graphically in Figure B.2. Since 
we are only concerned with the positions (x) v/ithin 


-P,< X <+p^ 

only (B12) and (B13) need be considered. These two equations and their 
graphic representation illustrate the fundamental parabolic nature of the 
acceleration limiting of any notion base simulator. 


I 

I 
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Superimposing the velocity limit on the parabolic limiting curves, 
the limiting envelope is as shown in Figure B.3. 




then X, 


= +v 


(Bio) 



I 


■Pi ^^desired^ 'Pi ^ ^’^current’^desired^^^ 

then “ ■ \/ ^^i^’^desired^Pi^ 

■Pi ^^desired*" -p^ & (x^urrent'^desired^^^ 

then x^jf^ij = -Vj^ 

^ 1 ^desired I ^ Pi ' ’^LIMIT ' ° 


(B19) 


(B20) 

(B21) 


The condition (x,,,,ent*^desired^ the sinful ator 

is moving towards the parabolic velocity limit boundary or away from it 
(equivalently, towards the center point or away from it). 


% 
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ure 5: Motion drive logic for the FSAA for longitudinal and pitch axis or lateral 
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Figure 6: Coordinate Systems [5] 



Centroid 

payload platform 



Figure 7: Vector Relationships for Actuator i [5] 
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Figure 9: Relation Between Vectors R and s 
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Flt’jre 10: Mutton Drive L 09 IC for Simulating an Sr.MS on the fSAA 
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